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First-principles  calculations  have  been  employed  to  study  clusters  of  Zr  embedded  in  Si02.  Stable 
complexes  are  found  with  four,  six,  and  seven  oxygens  coordinated  to  the  Zr  atom.  Consistent  with 
experiment,  the  higher  coordinated  complexes  are  the  most  stable.  These  also  have  a  higher  density, 
and  hence,  smaller  molar  volume.  This  smaller  molar  volume  provides  an  explanation  of  the 
increased  dielectric  constant  of  Zr^Si  ,  -x02  mixed-oxide  systems  for  small  amounts  of  Zr  (x 
<0.3).  An  unusual  sevenfold  coordinated  structure  is  described.  ©  2001  American  Institute  of 
Physics.  [DOI:  10.1063/1.1409590] 


The  technology  driven  search1  for  a  replacement  of  Si02 
in  applications  such  as  transistor  gate  insulators  and  memory 
cell  and  coupling  capacitor  dielectrics  has  lead  to  a  study  of 
a  wide  range  of  materials.  For  a  given  material  to  be  chosen 
it  should,  ideally,  have  a  large  permittivity  with  respect  to 
Si02  (say,  8  —  20-30),  have  a  large  band  gap  (to  avoid  elec¬ 
tron  or  hole  injection  from  the  gate  electrode  or  substrate), 
and  be  thermally  compatible  with  Si  (i.e.,  remain  amorphous 
and  avoid  a  reaction  resulting  in  the  creation  of  silicides  or 
roughened  Si/insulator  interfaces).  At  the  present  time  no 
clear  candidate  has  emerged  but  interesting  results  have  been 
obtained  with  Zr^Si^j-CK  (x<0.3)  mixed  oxides,2  4  where 
the  dielectric  constant  is  —11  for  x— 0.2.  The  dielectric  con¬ 
stant  increase  is  probably  too  small  with  respect  to  Si02  to 
merit  the  industrial  investment  necessary  to  change  technol¬ 
ogy  to  this  material.  However,  we  note  that  on  a  simple 
linear  extrapolation  basis  one  would  expect  the  dielectric 
constant  of  this  compound  to  be  —7  so  that  some  clear  and 
important  enhancement  is  occurring.  Note  that  enhanced  di¬ 
electric  constants  are  also  observed  in  other  mixed-oxide 
compounds5  such  as  Ta205/Ti02  so  that  this  effect  is  not 
limited  to  the  Zr02/Si02  system.  What  is  the  physical  origin 
of  this  enhancement  and  can  one  engineer  it  in  more  techno¬ 
logically  useful  materials? 

The  dielectric  constant  of  a  mixed-oxide  system  can  be 
treated  in  the  framework  of  the  oxide  additivity  rule6'7  where 
the  dielectric  constant  can  be  expressed  as 

e  =  {l+[8  7r/(3Vm)]  °tm{M  AxM'BpO Ay  vBq)}l  { 1 

—  [47t/(3  Vm)]am(MAxM'BP0Ay+Bq)},  (1) 

where  am(MA:M'BpOAy  ,  Hq)  is  the  molar  polarizability  for 
the  mixed-oxide  MxOy/M'pOq  (written  as  a  weighted  sum  of 
the  individual  molar  polarizabilities6'7)  and  Vm  is  the  molar 
volume  of  the  mixed-oxide  molecule.  Since  the  denominator 
of  Eq.  (1)  contains  a  negative  sign,  the  size  of  the  product  of 
the  molar  polarizability  and  of  the  molar  volume  with  respect 
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to  unity  are  clearly  crucial  in  determining  the  magnitude  of 
the  mixed-oxide  dielectric  constant.  In  the  following,  we  re¬ 
port  the  results  of  a  theoretical  investigation  of  the  behavior 
of  the  molar  volume  in  mixed-oxide  compounds  with  spe¬ 
cific  reference  to  the  Zr02/Si02  system.  We  demonstrate  that 
bonding  coordination  plays  a  crucial  role  and  that  in  fact 
high  coordination  leads  to  substantial  molar  volume  de¬ 
creases,  which  should  result  in  enhanced  dielectric  constants. 
The  importance  of  volume  effects,  as  explicited  in  Eq.  (1), 
was  neglected  in  an  earlier  explanation  of  the  origin  of  the 
enhanced  dielectric  constant.4  Furthermore,  in  the  model  out¬ 
lined  in  Ref.  4,  fourfold  coordination  was  extensively  evoked 
whereas,  as  we  shall  demonstrate,  our  conclusions  are  to  the 
contrary  that  higher  coordination  is  important  in  Zr02/Si02 
binary-oxide  compounds.  The  importance  of  all  these  factors 
will  become  evident  in  the  following. 

We  have  used  first-principles  theory  to  study  clusters 
simulating  a  Zr  atom  embedded  in  Si02.  Two  types  of  cal¬ 
culations  were  performed — density  functional  theory  (DFT) 
based  on  the  B3LYP  functional  using  the  North  West  Chem¬ 
istry  program8  and  second-order  many-body  perturbation 
theory  based  on  a  Hartree-Fock  reference  (MP2)  using  the 
Gaussian  atomic  and  molecular  electronic  structure  system 
program.9  To  allow  the  study  of  fairly  large  clusters,  the 
effective  core  potentials10  were  employed  in  all  calculations 
along  with  the  corresponding  valence  shell  basis  set.10 

Initial  clusters  considered  a  central  Zr  coordinated  to 
oxygen  neighbors  that  were  bonded  to  Si  next  neighbors.  The 
three  remaining  Si  bonds  on  each  Si  were  terminated  with 
hydrogen  atoms.  Complexes  with  different  numbers  of  oxy¬ 
gens  coordinated  to  the  Zr  were  found,  corresponding  to  sev¬ 
eral  different  local  bonding  schemes.  Stable  clusters  have 
been  identified  with  the  Zr  coordinated  to  four,  six,  and 
seven  oxygens.  The  simplest  clusters  are  the  four-  and  six- 
coordinated  systems  with  formulas  Zr(OSiH3)4  and 
Zr(0SiH3)4[0(SiH3)2]2 .  In  the  four-coordinated  cluster  the 
Zr  occupies  a  simple  substitutional  tetrahedral  Si  site.  Three 
different  six-coordinated  structures  have  been  found  which 
vary  only  in  the  relative  orientation  of  the  H3SiOSiH3 
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FIG.  2.  Stable  seven-coordinated  Zr  structure. 


FIG.  1.  Optimized  structures  with  Zr  coordinated  to  four  and  six  oxygens. 
The  central  atom  is  zirconium,  the  small  atoms  are  oxygens,  and  the  remain¬ 
ing,  large  atoms  are  silicons. 

groups.  Figure  1(a)  shows  the  most  stable  complex  for  four 
nearest-neighbor  oxygens  and  Fig.  1(b)  for  six  oxygens. 
Both  DFT  and  MP2  calculations  indicate  that  the  six- 
coordinated  structure  is  more  stable  than  the  four- 
coordinated  structure  plus  two  uncoordinated  H3SiOSiH3 
units.  Energetically,  the  difference  is  about  0.5  eV  in  the 
DFT  case  and  1.5  eV  in  the  MP2  case. 

From  examination  of  these  two  types  of  local  environ¬ 
ments  it  is  apparent  that  Zr’s  ability  to  coordinate  to  more 
than  four  oxygens  causes  a  significant  change  in  the  local 
network  [Figs.  1(a)  and  1(b)].  The  four-coordinated  system 
has  four  first-nearest-neighbor  Si  atoms  and  the  following 
layer  would  involve  12  oxygen  atoms.  The  six-coordinated 
system  has  eight  first-nearest-neighbor  Si  atoms  and  a  result¬ 
ing  24  oxygens  in  the  next-nearest-neighbor  layer — twice  the 
number  from  the  four-coordinated  system.  This  higher  local 
network  density  would  give  rise  to  a  smaller  molar  volume 
for  these  Zr  sites  and,  as  a  consequence,  give  larger  e  [Eq. 
(1)].  Based  on  the  above  models,  the  density  computed  in¬ 
side  a  sphere  with  radius  5.0  A  centered  on  the  Zr  atom  is  1.6 
times  larger  for  the  six-coordinated  structure  than  for  the 
four-coordinated  structure.  This  increase  in  density  for  high 
coordination  of  the  Zr  clearly  corresponds  to  a  decrease  in 
molar  volume  V m ,  and  therefore,  with  reference  to  Eq.  (1), 

an  increase  in  the  dielectric  constant  of  the  material. 

Downloaded  06  Feb  2007  to  129.238.237.96.  Redistribution  subject 


A  very  interesting  seven-coordinated  structure  has  also 
been  found.  Starting  the  calculation  from  a  structure  similar 
to  the  six-coordinated  one  with  a  formula 
Zr(0SiH3)4[0(SiH3)2]3  resulted  in  a  rearrangement  of  the 
ligands  to  produce  the  structure  shown  in  Fig.  2.  In  this 
optimized  structure  three  of  the  OSiH3  ligands  have  com¬ 
bined  with  the  three  H3SiOSiH3  ligands  to  form  three  un¬ 
usual  bidentate  ligands  (two  connections  to  the  Zr).  Each  of 
these  new  structures  contains  a  five-coordinate  Si  atom,  sug¬ 
gesting  the  presence  of  a  negative  charge  and  significant 
charge  separation  from  the  Zr  to  the  Si  atoms.  This  structure 
is  stable  in  both  DFT  and  MP2  calculations  but  not  in  the 
Hartree-Fock  method.  When  the  total  energy  is  compared  to 
the  six-coordinate  structure  plus  an  unattached  H3SiOSiH3, 
the  DFT  approach  indicated  the  six-coordinate  structure  is 
more  stable  but  the  MP2  calculations  favor  the  seven- 
coordinate  structure.  The  effects  of  embedding  this  cluster 
into  a  larger  network  will  need  to  be  included  to  accurately 
determine  its  relative  stability  and  further  calculations  are 
ongoing. 

The  further  increased  coordination  of  the  Zr  again  causes 
the  local  density  to  be  increased  (now  30  second-neighboring 
oxygens),  which  will  enhance  the  dielectric  constant  through 
a  decrease  in  V  m  .  The  density  inside  the  5.0  A  sphere  is  now 
1 .9  times  the  four-coordinated  density.  In  addition,  the  afore¬ 
mentioned  charge  delocalization  should  have  the  effect  of 
increasing  the  polarizability  term,  am ,  in  Eq.  (1) — also  in¬ 
creasing  the  dielectric  constant  for  this  type  of  system. 

It  is  worthwhile  at  this  point  to  compare  the  conclusions 
that  can  be  drawn  from  the  calculations  with  those  extracted 
from  experiments.  From  extended  x-ray  absorption  fine 
structure  measurements  on  xerogels"  it  has  been  ascertained 
that  the  O  coordination  around  the  Zr  atoms  is  3s  7  and  this 
for  concentrations  of  Zr01Sio.ci02.  For  trace  quantities 
(—2000  ppm)  of  Zr  in  a  variety  of  predominately  Si02-based 
glasses,  the  first-nearest-neighbor  O  coordination  has  been 
found12  to  be  6.  Finally,  over  a  range  of  concentrations  of 
ZrxSi  |  _r02  compounds,  the  O  coordination  has  again  been 

found13  to  be  6.  It  is,  therefore,  tempting  to  believe  that  the 
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zirconium  first-nearest-neighbor  O  coordination  in  the  mixed 
ZrOi  /Si02  oxides  is  6. 

The  calculations  presented  enable  us  to  conclude  that  the 
effective  molar  volume  is  directly,  and  rather  critically,  re¬ 
lated  to  coordination  of  the  Zr  atom  in  the  Si02  network. 
Therefore,  one  can  anticipate  that,  in  general,  dissolving 
high-coordination  atoms  into  low-coordination  networks 
such  as  SiOi  will  inevitably  lead  to  an  anomalously  rapid 
decrease  in  the  molar  volume/increase  in  the  molar  density 
around  the  dissolved  atom.  This  molar  volume  decrease  will 
manifest  itself  in  the  dielectric  constant  of  the  network.  At 
the  present  time  we  are  unable  to  absolutely  quantify  the 
density  changes  expected  in  the  “real”  Zr02/Si02  mixture. 
Further  calculations  on  the  Zr02/Si02  network  are  under¬ 
way  on  even  larger  clusters  to  (1)  study  the  stability  and 
energetics  of  forming  connections  in  the  a-Si02  materials 
and  (2)  to  add  additional  Zr  atoms,  i.e.,  extend  the  size  of  the 
clusters  to  include  more  than  one  Zr  atom. 
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